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Abstract 

A systematic inelastic neutron scattering study of the superexchange in- 
teraction in three different undoped monolayer cuprates (La2Cu04, Nd2CuC>4 
and Pr2Cu04) has been performed using conventional triple axis technique. 
We deduce the in-plane antiferromagnetic (AF) superexchange coupling J 
which actually presents no simple relation versus crystallographic parameters. 
The absolute spectral weight of the spin susceptibility has been obtained and 
it is found to be smaller than expected even when quantum corrections of the 

AF ground state are taken into account. 
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The copper spins properties of the insulating cuprates are of particular interest as they 
give insights into the microscopic description of the high-T^ superconductors. Undoped par- 
ent compounds of many high-Tc cuprates are usually described as Mott-Hubbard insulators. 
They exhibit an antiferromagnetic ordering below a Neel temperature ranging between 250 K 
and 420 K. This Neel state is well accounted for by a spin-| antiferromagnetic (AF) Heisen- 
berg model on a square lattice [Q. The following Hamiltonian, H = —JJ2<ij> SiSj where 
the sum is performed over spin pairs, is then used to describe the AF ground state where 
the most essential and generic parameter is the huge Cu-O-Cu superexchange interaction, 
J, within the CuO^ plane. J is usually determined by inelastic neutron scattering (INS) ex- 
periments which probe the dispersion relations of spin-wave excitations. The intraplane AF 
superexchange is then deduced from the measured spin-wave velocity c = 2SV2Z c Ja 
(where a is the square lattice constant, S=| and Z c ~ 1.18 represents quantum corrections 
of the AF ground state). Unfortunately, due to the large steepness of the in-plane spin wave 
dispersion (related to the large value of J ~ 100-150 meV), the spin- wave velocity is not 
easily deduced from INS experiments. Therefore, a precise knowledge of J is still needed in 
parent compounds of cuprates. Another essential magnetic parameter is the spectral weight 
of copper spin susceptibility which has been, so far, only reported in La2Cu04 The 
importance of these two parameters has been recently emphasized in doped materials as J 
is found to be renormalized compared to the undoped case and the spectral weight is shifted 
to lower energy [||,|J. 

Here, we present, by systematic neutron scattering measurements, the spin wave excita- 
tions of three different parent compounds of single-Cw02 layer cuprates. Especially, using an 
adapted focalised neutron scattering geometry, we are able to determine their spin-velocity 
with accuracy and to deduce J. Furthermore, we have determined the spectral weight of 
the spin susceptibility in absolute units and the perpendicular spin susceptibility, x±- X± 
can be also obtained as a consequence of sum-rules by applying the hydrodynamics relation, 
p s = (c/a) 2 x±_ P|JT5[], where p s is the spin-stiffness constant. We found that x± measured 



in neutron experiments is smaller than expected from this relation. This reduction of about 



30% is presumably due to covalent effects between copper d-orbitals and oxygen p-orbitals. 

High quality La 2 Cu0 4 , Nd 2 Cu0 4 and Pr 2 Cu0 4 single crystals of similar volume of about 
0.5 cm 3 were used. Neodymium and Prasedymium-based samples exhibit a Neel temperature 
around 250 K whereas the AF transition occurs just above room temperature, 320 K, in 
the Lanthanum-based sample ||. The samples were mounted with the reciprocal directions 
(110) and (001) within the scattering plane [these directions are referring to the tetragonal 
reciprocal lattice with Q = (h,h,q c ). We used the same axis in the case of orthorhombic 
La 2 Cu0 4 ]. Inelastic neutron scattering has been performed on the triple axis spectrometers 
IT and 4F1, installed respectively on thermal and cold source beams at the Orphee reactor, 
Saclay. The (002) reflection of Pyrolytic Graphite was used for both monochromator and 
analyser. No collimation was used and a filter (Graphite one on IT and Beryllium one on 
4F1) was placed on the scattered beam to remove higher order contaminations. 

A special scattering geometry || was used in order to align the resolution spectrometer 
ellipsoid along the AF line, i.e. the (001) direction. Namely, this focalisation allows us to 
separate counterpropagating spin-waves at relatively low energies as compared with standard 
geometries (7|||. We extend this technique down to 15 meV. For such a geometry, only one 
q c value is accessible for a fixed energy transfer and a fixed final neutron energy. To be 
powerful, this geometry also requires very good sample mosaicities. 

We now present q-scans (constant energy transfer scan) along the (110) direction in the 
three different monolayer cuprates: La 2 Cu0 4 , Nd 2 Cu0 4 and Pr 2 Cu0 4 . Figure [I] depicts q- 
scans measured at an energy transfer around 60 meV using the same experimental setup. The 
double peak structure is clearly seen in La 2 Cu0 4 and in Pr 2 Cu0 4 whereas only a flattened 
peak shape is observed in Nd 2 Cu0 4 . This difference emphasizes a larger spin velocity in 
Nd 2 Cu0 4 . In order to improve at low energy the determination of the spin velocity, we have 
applied in Pr 2 Cu0 4 this focalised geometry down to E — 14.5 meV, where a flattened peak 
shape is found (Fig. Q). Our data in Pr 2 Cu0 4 represent a clear improvement of a previous 
measurement JTG[] . 

Here, we focus on the low energy part of the spin wave spectrum in the limit where the 
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dispersion relation for AF magnons is linear (hu << 2Z C J). However, at low energy, the 
magnon spectrum exhibits gaps which are either related to planar anisotropy or to interlayer 
interactions J?]. The usual linear relation is thus only recovered for energies slightly larger 
than these gaps. Due to the large intraplane superexchange interaction in cuprates, this 
condition is fulfilled for energy above ~ 12 meV (see Fig. |2|). Above this energy, the spin- 
wave neutron cross section per formula unit can be written in terms of the dynamical spin 
susceptibility fll|Jl2"l| , x(Q)^)> as 



d 2 a _ 2 F 2 (Q) 1 (3 Q 2 C) Im X (Q,uj) ^ 



dVtduj u vr(^ B ) 2 2 v Q 2J l -exp(-fuv/k B T) 

where Tq=0.292 barns, -F(Q) is the atomic form factor of Cu 2+ spins g ~ 2 is the Lande 
factor for copper spins, and Q c = ^-q c is the component along the (001) direction of the 
scattered wavevector, Q. For an AF single layer cuprate, the imaginary part of dynamical 
susceptibility of the low energy spin wave excitations is given in absolute units by |12| 



, 2 V2 t 



Imx(Q,u) = S7cZ x Z c (gfi B ) —5[u-cq] (2) 

where q is the in-plane wavevector component along the (110) direction referred to the AF 
wavevector. The quantum corrections associated to the perpendicular susceptibility Z x , 
is included. The convolution product of the Gaussian resolution ellipsoid by the spin-wave 
cross section (|I|) with the spin susceptibility (§) gives i) the dispersion relation of magnons ii) 
the spectral weight of Imx- The q-scans have been fitted by this convolution product with 
4 fitting parameters: the magnon in-plane wavector q, the amplitude of Imx an d a sloping 
background. We note that the observed experimental q- width along the (110) direction 
merely corresponds to that of the resolution. 

In Pr 2 CuO"4, the in-plane magnon dispersion is reported in Fig. Q over a wide energy 
range. As expected, a linear dispersion typical of AF excitations is found with a slope which 
is the spin wave velocity, c= 0.80 eV.A. Comparison of the different q-scans (fig. |I]) gives 
0.85 eV.A for La2Cu04 in agreement with a previous determination by high energy neutron 
experiments [§] and c=1.02 eV.A for Nd2CuO"4 (see Table (0])). The magnon wavevector, 



and so the spin velocity and the AF intraplane super exchange, are then found larger for 
Nd 2 Cu0 4 by about 20% as compared with the two other systems. 

The spin susceptibility in absolute units has been experimentally estimated by a stan- 
dard calibration [|J using acoustic phonons, whose dynamical structure factor is known 
by lattice dynamics. The magnetic part has been measured from high energy scans (Fig. 
[I]) as well as non-resolved low energy q-scans. In order to compare the observed spin 
susceptibility in absolute units with its theoretical predictions [Q, we calculate the aver- 
age of (@) over the two dimensional (2D) q-space perpendicular to the (001) direction, 
X.2D = I rfq2D-fw\;(Q) ^0/ / dq.2T>- In our experimental energy range, \2D is almost indepen- 
dent of energy: X2D — S(gnB) 2 Z x /2J . Values for \2D are listed in Table (|). In La 2 Cu0 4 , it 
compares well with two previous measurements |2|||. On the one hand, Itoh et al. have 
reported an effective value of S=0.17 which is reduced from the spin number, S=l/2. That 
agrees with our observed spin susceptibility, 2.7 p 2 B /eV (see Table (|)), which is reduced 
by the same factor from the classical spin susceptibility (without quantum corrections), 
X2D S — S(g^B) 2 /2J = 7.5 /ig/eV. On the other hand, Hayden et al. || have obtained 
X.2D = 2.5 /z^/eV which agrees in errors with our value [TJ[. 



The perpendicular susceptibility, x±> deduced from our INS measurements is then ob- 
tained by applying the relation x± — X2d/4S , ((? / ub) 2 |l| and listed in Table (|). x± can be 
independently deduced from the spin stiffness, p s , applying standard hydrodynamics rela- 
tion in the Heisenberg model (see Table (|)). Let us recall that the spin-stiffness constant 
has been obtained in the Heisenberg model from the two-dimensional correlation length £ 2 d 
above the Neel temperature as, £ 2 d oc exp(|^) Hl5| , being itself measured using energy 
integrated neutron scattering [|17l , |i"8f . Surprisingly, the value of x± deduced from p s is found 
to be systematically smaller than the one measured in INS experiments even when quantum 
corrections of the AF ground state are considered. This discrepancy of about 32% for the 
three compounds is likely due to the covalence of copper d-orbitals with oxygen p-orbitals 
TJ|. Reducing the absolute scale of the atomic form factor, such effects can explain the 



diminution of the inelastic spectral weight of the spin susceptibility as well as the low tern- 



perature ordered magnetization value 0. Consequently, the spectral weight of Imx does 
not solely determine the quantum corrections for the spin susceptibility 

We now deduce J as well as the quantum corrections. Since there are more unknown 
parameters that the measured ones, we need to use theoretical estimation for one parameter. 
Among the measured magnetic parameters, the spin wave dispersion curve is presumably 
the less altered by frustration effect and disorder ||20|| . The quantum correction to the spin 
wave velocity Z c estimated from different theoretical approaches fl]]2(| likely converges to 



a best value of Z c = 1.18 HTfJ. J is then confidently deduced from the spin wave velocity 
using this value (see Table (|)). Two other parameters are related to J. On the one hand, 
the spin-stiffness constant is usually modelled as p s = Z Ps JS 2 [^0| (where Z Ps accounts 
for quantum effects). On the other hand, a high frequency broad peak is observed in Ra- 
man scattering which is likely interpreted as two-magnons processes with opposite momenta 
pl| , p2| . By means of series expansions technique |[23|| , the moments of the Raman intensity 
(the frequency of the spectrum maxima Mi as well as lineshapes) have been related to J, for 
instance M\j J = 3.58. The quantum corrections for the spin stiffness Z Ps , the perpendicular 
susceptibility Z x , and the ratio between the first Raman scattering moment and J have been 
obtained and also listed in Table (|). 

Surprisingly, only the quantum corrections found in La2Cu04 are in agreement with 
the theoretical predictions IJ either based on series expansions [0,^] or based on 1/S 



expansion linear spin-wave theory fl6| : Z Ps = 0.72 and Z x = 0.51 and ojr/J=3.58. The 



two other systems display larger quantum corrections for p s and x± ma y be related to their 
different low energy spin excitations 0. An even larger discrepancy is observed for the 
spin pair Raman scattering measurements. Consequently, the neutron measurements which 
determine p s as well as the light scattering experiments only give a rough estimation of J. 

We now relate the copper spin intraplane superexchange determined by INS with the 
crystallographic distances between copper atoms (Figure 0). Clearly, J does not exhibit a 
monotonous dependence versus the bonding Cu-O-Cu length in contrast to what could be 
expected. This outlines that the classical superexchange theory being only related to the 



Cu-O-Cu bonding is a too simple description. Moreover, it has been recently stressed that 
the large enhancement of J is actually caused by another structural unit, namely the Cu- 
0-0 triangle |p3| . Empirically, one can distinguish distorted tetragonal lattice and perfect 
square one. Indeed, J appears to decrease sharply with the distances between copper atoms 
in Nd 2 Cu0 4 and in Pr 2 Cu0 4 (both having the T'-phase, i.e. linear Cu-O-Cu bonding). Note 
that the largest J is found in Nd2Cu04 where the Cu-0 distance exactly corresponds to the 
sum of copper and oxygen ionic radius. The two other systems do not belong to the same 
family as the bonding Cu-O-Cu is not linear: it is distorted perpendicular to the plane in 
YBa 2 Cu306+a; 0, or even in both directions in La2Cu04 |25[ due to the tilt of the CuOq 
octahedra. Therefore, J turns out to be extremely sensitive function of Cu-O-Cu bonding 
angle. 

In conclusion, by means of inelastic neutron scattering experiments using conventional 
triple-axis technique, we deduce J and the quantum corrections of the AF ground state 
in undoped monolayer cuprates. The in-plane antiferromagnetic superexchange coupling 
J does not exhibit a monotonous behaviour versus the bonding Cu-O-Cu length. The 
absolute spectral weight of the spin susceptibility is smaller than expected from quantum 
corrections JTJ, likely due to covalent effects. These results provide a necessary ground for 
the understanding of antiferromagnetism in the high-Tc superconductors. 
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TABLES 



Parameter 


T N 


c 


X2D =>" 


XX (INS) 


2irp s => 


X± 


J 




Z X — Z Ps/ Z l 




Units 


K 


meVA 




eV" 1 


meV 


eV" 1 


meV 








Errors 




± 20 


± 0.4 


± 0.05 


± 5 


± 0.04 


± 3 


± 0.05 


± 0.04 




La2CuC>4 


320 


850 


2.7 


0.34 


150 a 


0.48 


133 


0.72 


0.52 


3.5 C 


Nd 2 Cu0 4 


246 


1020 


1.8 


0.22 


137 6 


0.33 


155 


0.64 


0.46 


2.5 C 


Pr 2 Cu0 4 


252 


800 


2.3 


0.29 


114 6 


0.44 


121 


0.6 


0.43 


3.1 d 



TABLE I. Magnetic parameters in three undoped single layer cuprates. The value of the spin 
stiffness has been deduced from previous energy-integrated neutron scattering experiments: a from 
[17], b from [18,10|. ojr is the first moment of the Raman scattering data: c from d from [22]. 
Note that T/v is not simply related to J due to the 2D character of the magnetic interactions in 
cuprates M. 
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FIG. 1. q-scans across the magnetic line around %uj ~ 60 meV in three different monolayer 
undoped cuprates. Typical counting time is 1 hour per point. Full lines correspond to the convo- 
lution product of the Gaussian resolution ellipsoid by the spin-wave cross section (D with the spin 
susceptibility (pi). 
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FIG. 2. Left: q-scan across the magnetic line at Tluj = 14.5 meV in Pr2Cu04 (see Fig. I 

for details). Right: In-plane magnon dispersion in Pr2Cu04. At low energy, the degeneracy 

between out-of-plane and in-plane spin components is removed due to planar anisotropy leading 

to an out-of-plane spin gap of about 8 meV Q. Above ~ 12 meV, both spin components become 

very rapidly indistinguishable with increasing the energy. Open circles correspond to a previous 

measurement |1Q| . 
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FIG. 3. In-plane superexchange interaction versus Cu-O-Cu bonding length in different 



cuprates. The value for the bilayer system YBCO is from 
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